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\yr{ ' Abstract 

O ! 

<^ , We calculate the cosmic microwave background (CMB) anisotropy 

CN | spectrum in models with millicharged particles of electric charge q ~ 

O ■ 10~ 6 — 10 _1 in units of electron charge. We find that a large region of 

the parameter space for the millicharged particles exists where their 
effect on the CMB spectrum is similar to the effect of baryons. Using 
■ WMAP data on the CMB anisotropy and assuming Big Bang nucle- 

. osynthesis value for the baryon abundance we find that only a small 

J> | fraction of cold dark matter, O mcp /iQ < 0.007 (at 95% CL), may con- 

sists of millicharged particles with the parameters (charge and mass) 
from this region. This bound significantly narrows the allowed range 
of the parameters of millicharged particles. In models without para- 
photon millicharged particles are now excluded as a dark matter can- 
didate. We also speculate that recent observation of 511 keV 7-rays 
from the Galactic bulge may be an indication that a (small) fraction 
Q_i| of CDM is comprised of the millicharged particles. 

QJ , Search for particles carrying small but non-vanishing electric charge (mil- 

licharged particles) has long history. If observed, millicharged particles would 
either cause serious doubts on the concept of Grand Unification or imply the 
^ ■ existence of a new massless gauge boson — paraphoton 012] • Furthermore, 

the existence of millicharged particles would hint towards processes with 
apparent electric charge non-conservation, like electron or proton decay to 
"nothing" 0. 

There are various constraints on the parameters (charge and mass) of 
millicharged particles, coming from collider and laboratory experiments and 
from cosmology and astrophysics (see, e.g., Refs. [11001111313 for the latest 
results and Refs. jTUl HI] for reviews), see Fig. [TJ Interestingly, reported 
bounds did not exclude a possibility [12] that a significant part (or even all) 
of the cold dark matter (CDM) is comprised of millicharged particles. 
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The constraints on the parameters of millicharged particles are somewhat 
different in theories with and without paraphoton. Without paraphoton, 
two domains in the parameter space of millicharged particles are allowed. 
The first one corresponds to heavy particles with tiny electric charge (left 
upper corner in Fig. [TJ which would never be produced thermally in the 
early Universe. This region is far beyond the reach of future collider and 
laboratory experiments. In the current Letter we are concerned with another 
region. This is a narrow window of relatively light particles with masses 
M ~ 10~ 3 10 2 GeV and charges q ~ 10 _6 -i-10 _1 . Larger charges q > 0.1 are 
ruled out by limits on the cosmic ray fluxes of fractionally charged particles 
[llj. It is worth noting that particles with q > 0.2 are also excluded by 
measurements of the width of Z-boson (cf. Ref. [3]) if one makes use of the 
latest data [T3] . 

In model with paraphoton, and for not very small values of the parapho- 
ton coupling constant a', millicharged particles annihilate mainly into pairs 
of paraphotons. As a result, their annihilation in the early Universe is more 
efficient and the cosmological bound coming from the relic abundance de- 
pends on the value of a' and generically is less restrictive (see Fig. ^) than 
in the model without paraphoton. 

It was noted in Ref. ^1] , that there is a part of the parameter space for 
the millicharged particles where they do not decouple from the acoustic oscil- 
lations of the baryon-photon plasma at recombination, and it was suggested 
that the effect of these particles on the CMB anisotropy spectrum may be 
similar to the effect of baryons. The purpose of this Letter is, using the recent 
precise CMB data from WMAP ^2], to set an upper limit on the fraction 
of millicharged particles in CDM and to narrow the allowed window for mil- 
licharged particles. Assuming the standard Big Bang nucleosynthesis (BBN) 
value for the baryon abundance, Vt^h^ = 0.0214 ± 0.0020 jTB], we arrive at 
the following constraint on the millicharged particle abundance, 

tt mcp h 2 < 0.007 (95% CL) , (1) 

if millicharged particles are coupled to baryons at recombination. The latter 
condition is satisfied on the right of the dark solid line in Fig. ^ We see 
that the upper limit ([!} apllies in the whole allowed window for millicharged 
particles in models without paraphoton. It is also worth stressing that in 
models with paraphoton the domain of applicability of the upper bound 
does not depend on the value of a'. Using the Lee- Weinberg formula [17] 
for the relic abundance one then translates the upper bound ((TJ) into the 
lower limit on the annihilation cross section of millicharged particles. As a 
result, the upper bound (JTJ) excludes most part of the allowed window for the 
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Figure 1: The exclusion plot in the parameter space for millicharged par- 
ticles. Light grey area is excluded by accelerator experiments and BBN. 
Dashed region is excluded by the relic abundance of millicharged particles in 
models without paraphoton. Part of this region above the dash-dotted line 
is excluded by the relic abundance in models with paraphoton (assuming 
a 1 = 0.1). On the left of the dotted line millicharged particles cannot be 
thermally produced in the Early Universe. Dark grey and dashed dark grey 
areas are the previously allowed regions which are now excluded by Eq. 
in models without and with paraphoton, respectively. On the right of the 
dark solid line millicharged particles are coupled to baryons (see Eq. (JHJ) in 
the text). 

millicharged particles in models without paraphoton, leaving a small allowed 
region with masses in the range m ~ 10 _1 — 10 GeV and charges in the 
range q ~ 10~ 3 — 10 _1 . In models with paraphoton, the bound translates 
into meaningful limit on the parameters of millicharged particles as shown 
in Fig. [U the excluded region depends on the value of a'. 

Let us proceed to the derivation of the upper bound ((T}. To calculate the 
CMB anisotropy spectrum we adapt the CMBFAST code [IB] which solves 
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numerically the set of kinetic equations ^H] for the linear perturbations in 
the primordial plasma. To take into account the presence of millicharged par- 
ticles we extend this set by adding the kinetic equations for the millicharged 
component, modify the equations for the baryon component to take into ac- 
count the elastic scattering off millicharged particles and include millicharged 
component contribution to the energy-momentum tensor. The rest of the 
perturbation equations are the same as in Ref . • The Compton scattering 
off millicharged particles is negligible, since the corresponding cross section 
is suppressed by the fourth power of the charge q. 

We work in synchronous gauge and consider primordial plasma in the 
expanding Universe with scale factor a(r) (where r is a conformal time) 
normalized to unity at present time. Let Tf, pf, Vf be the temperature, 
density and velocity of the f-th component of the plasma. In particular, 
f = e,b, 7, mcp for electrons, baryons, photons and millicharged particles, 
respectively. In what follows bar denotes space averaging. The standard 
variables describing fluid perturbations are Sf(k, r) = \ff(k, r) — P/(t)]/ P/(t) 
and 8f(k,r) = ikjvj(k,r) where ki is conformal momentum. 

Before recombination, the interaction between non-relativistic electrons 
and protons is strong enough to ensure that electron and baryon components 
have equal velocities, 9 e = 6b- This makes it possible to use tight coupling 
approximation and consider electrons and protons as single baryon fluid [23] . 
Then the set of equations for baryons and millicharged particles reads (cf. 
Ref. H2|) 

8b = -Ob - \K 

0\ = --0 b + c 2 s k 2 S b + ^an e <7 T (0 7 - 6 b ) 
a 6pb 

-\- ClT mcp^"mcpi,0mcp Ob)-, (2) 



u mcp u mcp 2 

f) — --B +r 2 k 2 r) 

u mcp u mcp ~ ^s^mcp u mcp 

+ ar mC pf2;,(06 mcp ), (3) 

where h is the longitudinal metric perturbation, dot stands for derivative 
with respect to the conformal time r; c s , c Stmcp are the sound velocities in 
the baryon and millicharged components, n e is the number density of elec- 
trons and T mcp is the velocity transfer rate for millicharged particles due to 
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scattering off baryons and electrons. The latter is given by 



x=e,p ' 



M,x 



(4) 

Vm,x=0 



where brackets stand for thermal averaging, VM, e ( P ) is relative velocity and 
Avm is velocity transfer in a single process of scattering; dou.e (g?cm, p ) is 
the Rutherford cross section for millicharged particles scattering off electron 
(proton). 

The Rutherford cross section is singular at zero scattering angle, but due 
to Debye screening the integral in Eq. (j3J) is cut at the value of the scattering 
angle equal to the Debye angle, do — \j2itan e /T 2 m e . As a result one arrives 
at the following expression for the velocity transfer rate in the case of thermal 
equilibrium 



Ay2Tia 2 q 2 p crit . .. . , . 

TmCP = 3Mm p aV 2 (r)T^ ' * ^ ' ( ^ + ^ ' (5) 

where fiM,e( P ) is the reduced mass of a millicharged particle and electron (pro- 
ton), a is the fine structure constant and To ~ 2.726 K is the present CMB 
temperature. It is straightforward to generalize Eq. (JHJ) to nonequilibrium 
case when electrons, protons and millicharged particles have different tem- 
peratures. In that case the value of T mcp is larger than the one given by 
Eq. (j3J), so Eq. (jSJ) may be used as a lower estimate of the interaction rate, 
which is sufficient for our purposes. 

We solve the system of the kinetic equations starting from the early mo- 
ment of time Tj and using inflationary initial conditions jTH] , 

3 1 

&mcp — $b = J&y = ~~^p{kT,j) , 
1 .3 



Qracp — $7 — Qb — ~~T^Ck(kTj 
lo 



where constant C determines the overall normalization. 

The last terms in the r.h.s. of Eqs. (j2J), (jSJ) tend to equalize the velocities 
of the baryon and millicharged components of the fluid. This results in the 
kinetic equilibrium, 9 mcp = 9b, provided T mcp is large enough. In this case, the 
perturbation equations are difficult to solve numerically, because the kinetic 
relaxation rate for the interaction of millicharged particles and baryons is 
much larger than the rates of other processes. To deal with this situation 
we make use of the zeroth order tight coupling approximation adopting the 
method of Peebles and Yu Namely, we expand Eqs. O), © to the zeroth 
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order in , setting 9 b = 9 mcp = 9. Then we exclude T mcp from Eqs. (J2J), 
and arrive at the following equation 

e = --e+ o "o ^°^(g 7 - e) 

The CMB spectrum obtained in this approximation agrees with the solution 
of the original set of equations at the level of one percent provided that 

^mcp{'Trec)(y^b ~t~ ^mcp) H (T rec ) ^ 250 , (7) 

where r rec is conformal time at recombination and H(T rec ) is the Hubble 
parameter. In what follows we discuss a region of parameters where tight 
coupling condition (JJJ) is satisfied (in particular, this region covers the whole 
allowed window in models without paraphoton) and comment on the rest of 
the parameter space in due course. 

To compare the results of our simulation with the CMB data we consider 
the flat ACDM model with the number of massless neutrino species N v = 3. 
We perform a scan over the space of the cosmological models by varying 
parameters from minimal to maximal values as given in Table d All priors 
are at 95% CL. We assume helium fraction at the moment of recombination 



parameter 


rain, value 


max. value 


step 


reference 


&CDM 


0.2 


0.4 


0.01 


PDG 


13 


h 


0.64 


0.79 


0.01 


PDG 




n s 


0.8 


1.2 


0.01 






0.0194 


0.0234 


0.0005 


BBN [HI] 


^"mcp 





0.020 


0.001 





Table 1: The ranges of the cosmological parameters used in simulations; we 
never have a good fit outside the region n s 6 [0.8, 1.2], and for Q mcp > 0.02. 

Yjj e = 0.24. We have checked that reionization effects are irrelevant here, the 
reason being that reionization affects the spectrum only at the lowest values 
of the multipole moments. 

For each cosmological model we calculated likelihood to the WMAP 
data As a result we arrived at the limit given by Eq. (0), which means 
that no model exists in the considered parameter range with larger values 
of Qmcpho and likelihood better than 5%. We also checked that additional 
CMB data from CBI |2U and ACBAR [22] do not improve the limit JTJ). 
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Figure 2: Distribution of models in the (f^/iQ, ^mcp^o) pl ane - Crosses and 
dots denote models agreeing with data at the la and la CL, respectively. 
The bold line illustrates the degeneracy of the CMB anisotropy spectrum. 
Two dotted lines show the range of Qbh% allowed by BBN. 



Using the Lee- Weinberg formula we translate the limit given by Eq. 
into the bound on the parameter space. The corresponding excluded areas 
are shown in Fig. ^ In models with paraphoton we extend our analysis to 
the region where the tight coupling condition is not satisfied. We checked 
that actually the upper bound (JTJ) applies provided that the following less 
restrictive condition holds 

r mcp (x rec )(fi b + n rncp )H(r rec y l > 2.5 . (8) 

This inequality is true on the right of the black solid line in Fig. ^ 

Coming back to the tight coupling regime, let us illustrate the upper 
bound pictorially by generating a set of models with random parameters 
in tight coupling regime, assuming uniform distribution of models in the 
parameter space within ranges given in Tabled For each model we calculated 
likelihood to the observed CMB spectrum and plot the resulting distribution 
of models in the (f^/ip, fi mcp /ig) plane in Fig- El One observes that the CMB 
spectrum is approximately degenerate along the straight line 

Q b h 2 = 0.022 - 1.1 Q m cph 2 . (9) 

This degeneracy is in agreement with the expectation of Ref. [Hj, that the 
effect of millicharged particles is similar to that of baryons. Thus in models 
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Figure 3: Two different CMB anisotropy spectra compared with extended 
WMAP dataset. Solid line represents the best fit model without millicharged 
particles, VL b h^ = 0.022. Dashed line corresponds to model with VL h h^ = 
0.0U,n mcp h 2 = 0.007. 



with millicharged particles, the CMB data determine actually the sum 
(n b + n mcp )h 2 = 0.022±0.001 (68% CL). Combining this value with the lower 
limit Qbh>Q > 0.019 from BBN one arrives at the upper bound very similar to 
Eq. (P). This serves as a qualitative explanation of our result. 

Another illustration of the approximate degeneracy is shown in Fig. El 
where two CMB anisotropy spectra calculated for different models on the de- 
generacy line are shown. One observes that the two spectra almost coincide 
in the region of the first and second acoustic peaks. However, the degener- 
acy is no longer present at higher multipoles. This is due to the fact that 
the electroneutrality of the plasma implies that the electron number density 
is proportional to the baryon density. Hence, replacing a certain amount 
of baryons by millicharged particles results in the enhancement of the Silk 
damping at small scales. With future precise data for high values of I, one 
will be able to set a constraint on the value of VL mcp using the CMB data 
only, without reference to BBN results. To check this we created a simu- 
lated dataset, which contains the same values of I as in the WMAP data 
up to I = 500, and then with the step Al = 50 up to I = 1600. The 
CMB anisotropy spectral coefficients Cis were taken from the best fit ^Hj 
to WMAP data. The error bars for these coefficients were assumed to be 
equal to cosmic variance. Repeating the above procedure for this dataset we 
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obtained that an upper limit Q mcp hQ < 0.003 can be placed. Further improve- 
ment of this limit turns out to be impossible due to the new approximate 
degeneracy, 

Q b h 2 = 0.022 - 0.65 fl mcp h 2 Q , 

n s = 0.94 + 8.0 Q mcp , (10) 

arising at smaller values of Q mcp hl. 

To conclude, we note that when translated into the parameter space, 
the limit (JTJ) is especially interesting for the models without paraphoton, 
where it excludes most part of the window with not very heavy particles and 
substantial electric charges. To completely close the window, sensitivity to 
millicharged particle abundance at the level of flmcphy ~ 3 ■ 10~ 4 would be 
required, which cannot be achieved with future CMB experiments due to 
the degeneracy (jlOJ) . Determination of the baryon abundance from the BBN 
is not accurate enough to improve the situation. Hopefully, the rest of the 
window will be explored by future accelerator and/or laboratory experiments. 

Finally, recently it was suggested 23 j that the flux of the 511 keV 7-rays 
from the Galactic bulge detected by the INTEGRAL satellite j21] may be 
explained by the annihilation of the ~ 1 4- 100 MeV dark matter particles 
into e + e~ pairs, provided their annihilation cross section af3 and abundance fl 
satisfy (jf) (^ff^) 2 - lO"^ 4 - 5 ). Intriguingly, this condition holds 

in the left corner of the parameter space for millicharged particles without 
paraphoton allowed by Eq. (0) (say, for q = 3 ■ 10~ 3 , M = 100 MeV). 
One is tempted to speculate that the observation of the 511 keV line is an 
indication that a (small) fraction of CDM is comprised of the millicharged 
particles. This possibility will be checked by the future CMB data. 
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